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SECTION 1 


INTRODUCTION AND SUMMARY 

This Is the final report under Contract NASW-2222 sponsored by the 
Planetary Geology Program at NASA Headquarter (Code SL) . The objective of 
this program was to investigate and illustrate the application of radiometric 
analyses in the Interpretation of Mariner VI and VII imagery. We felt that 
the following t)q)es of problems could be addressed using this Imagery; 

(1) Evaluation of local reflectivity changes, 

(2) Augmentation of geologic mapping, 

(3) Discrimination of atmospheric phenomena, and 

(4) Investigation of polar cap structure. 

The use of conventional radiometric techniques requires the ability 
to Convert sensor response into surface radiance, that is radiometric calibration 
of the Mariner TV systems. In order to Interpret the surface radiance, 
knowledge of the bi-directional reflectivity function and the atmospheric 
effects (if any) are required. To circumvent the requirement that these 
functions, i.e., the radiometric response and bi-directional reflectivity, 
be known with a high degree of accuracy an Interpretation technique based upon 
reflectivity ratios was employed during this study. Alternate near encounter 
frames from the wide angle Mariner cameras were taken through different 
spectral filters (red, green or blue) . Consequently, the overlap between 
successive frames show the same area of the martian surface in two spectral 
bands. We call these areas the bi-band coverage areas. 
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A comprehensive review of previous and subsequent efforts to the 
present study was made and is included in this report. In performing this 
survey it was obvious to us that the terminology and definitions employed in 
discussing radiometric investigations (both inter-planetary and earth-based) 
is quite divergent and confusing to the reader. This lack of precision makes 
the comparison of results difficult. Consequently the next section of this 
report proposes a radiometric model which Incorporates the definitions and 
terms used in both astronomy and electro-optics. 

As noted above, the interpretation of radiometric data can require 
knowledge of the bi-directional reflectivity or photometric function. Con- 
sequently in Section 2 we have reviewed the definition of this function and 
discuss several parametric models for its representation. We also have 
included the parametric values obtained by several authors in fitting observed 
radiometric data to the functional representations for Mars. 

Section 3 reviews some of the radiometric analysis techniques 
employed in geology. The development of topographic information including 
crater distribution functions and topographic profiles, the use of global 
reflectivity in defining geologic units and in establishing surface composition 
and the Interpretation of local reflectivity differences are discussed. 

Our analysis on Mariner VI and VII imagery to interpret local 
variations in reflectivity by using reflectivity ratios is presented in 
Section 4. We identify the areas where bi-band spectral data is available 
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and present specific data in support of the interpretation methodology used. 
The results of our analysis show that light and dark markings interior to 
crater floors observed In Mariner frames 6N11 and 6N13 acquired through 
green and red filters respectively are caused by differences between particle 
sizes between the dark and light areas. Based upon laboratory spectral 
reflectivity data the most likely size groups for the particles in the two 
areas is consistent with that proposed by others. Global reflectivity 
differences between dark and light areas obtained by other authors support 
larger particle size differences. Based upon our data and that provided by 
other authors, we concluded that the most likely composition of this local 
area is limonlte stained p)nroxene, the latter being a basic silicate, common 

j. 

in meteorites . Desplte'the evidence for basaltic flows on Mars, olivine 

i ' . i , ' 

which is associated with such flows was ruled out as a major constituent in 
this area based upon its relatively high normal reflectivity. 

1 ' . 

i 

The application of reflectivity ratio analysis to geologic mapping 

i ■ - f ’ 

i and the study of atmospheric phenomena was Inhibited by the limited amount of 

I I 

data available. It was anticipated that the Mariner IX mission would provide 

! ; ■ 

a wealth of bl-band coverage areas on the martian surface and analysis of this 

data was anticipated. Unfortunately the filter wheel became inoperative on 

i 

1 

Revolution 118 subsequent to the subsiding of the global martian dust storm. 

As a consequence, no significant bi-band coverage area of the martian surface 
was obtained during the Mariner IX mission and this portion of the study was 
deleted. 
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SECTION 2 


THE RADIOMETRIC MODEL 

While the primary objective. of this study was not the development 

of basic radiometric (or photometric) principles or techniques but their 

application to planetology, an understanding of these principles had to be 

achieved in order to assess the efforts of related research. Unfortunately, 

the historical development of radiometric science in astronomy and in 

electro-optics has not had complete commonality- -leading to a diversity of 

definitions and structure that can be confusing to the engineer not versed 

in both disciplines. This can result in the misinterpretation of the data 

presented in other research efforts, Nlcodemus (1967) has noted: 

"Radlometry and particularly its overshadowing subdivision, 
photometry, are embarrassed by diversity in nomenclature. 

Careful attention to the definitions of all radiometric 
terms , symbols , and units both by authors and readers is 
needed to avoid confusion and misunderstanding." 

In this section of our report we present a proposed radiometric model and 

provide definitions that include concepts from both astronomical and 

i 

electro-optical radlometry as a basis for standardization in this and 
future efforts in terrestial geology as well as planetology. 

An Important distinction is made depending on the spectral response 
of the detector in radlometry. Terms such as luminance, illuminance and 
brightness are used in photometry where the detector has a spectral passband 
equal to that of the eye (photopic response). Although this is quite 
precise the prefix "photo" is frequently used when this is not the case. 
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An example In astronomy and planetology Is the definition of "photometric 
function" which can be measured over a narrow visible or IR spectral band 
and is rarely measured with a photopic response. The model and terminology 
presented below is based upon a summary of radiometry (Nicodemus, 1967). 

1 

The geometry of the incident and reflected light from a planar 
surface (,in the x-y plane) is shown in Figure 1. Reflectance of an opaque 
surface is a function of the Influx direction through angles B; and if): 
and the efflux direction through angles and . The basic quantity 


Z 



X 


Figure 1. RADIOMETRIC MODEL GEOMETRY 
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is the spectral radiance , emanating from the surface and is 

expressed in units watts/m *sr*^m. The wavelength dependence is removed 
by integrating over the bandwidth of the incident light and the spectral 
response of the detector. The explicit dependence on wavelength is 
suppressed below with the understanding that the terms imply a particular 
bandpass . 


The reflectance (dimensionless) of the surface is defined as the 
ratio of the reflected radiant power to the incident radiant power 

_ Pr_ /Nr(6rA)^-0-'^ ^ 

Pi /tJiie.jjjji'. f^i(&i,<t>i) 

where the projected solid angl S iv\Q coi6ei0j(P. Both the incident 

radiance and the reflected radiance are functions of their respective 
direction angles (0,(P) . The reflectance of the surface as defined above can 
change as the receiver geometry is changed. Clearly we do not generally want 
to adopt such a definition. An exception occurs when the detector measures 
all of the reflected radiation, l.e., that reflected into the hemisphere, in 
which case we obtain the diffuse reflectivity or the Bond albedo discussed 
below.* 


The reflecting properties of the surface are more appropriately 

described by the bidirectional reflectance function defined by 

*The term reflectivity implies decimal fraction while reflectance is expressed 
in percent. For historical interest we note that the Bond albedo was intro- 
duced in 1861 . 
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( 2 ) 


pYS; <!>■ e, dr) = (Sr, dr) 

where Hi (6i,d is the incident irradlance in watts/m^. Note that ^ has 
dimensions of steradlans Combining equations (1) and (2) the diffuse 
reflectivity (or Bond albedo) is given by: 


ej 


- -feXeiAM) Jn. 


( 3 ) 


where "h" signifies integration over the hemisphere. For a perfectly diffuse 
(Lambertian) surface p ^Is independent of receiver coordinates and 



(4) 


The normal reflectivity (or normal albedo ) , P , is the value of the 
bidirectional reflectance function when 0 ; =0 = 0 and the Influx and efflux 
directions are both normal to the surface; 

p = (O.O) 

Hi (0,0) 

In practice the normal albedo can be computed by measuring the radiance of 
the center of the planet (or the sub-earth point) near opposition and 
calculating Hi • 


If we normalize the bidirectional reflectance to 1.0 at zero influx 
and efflux zenith angles Equation (2) becomes: 
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( 6 ) 


where ^ Is the normalized function. In planetology and astronomy ^ is 
referred to as the photometric function even though it can be measured over 
non-photoplc spectral bandpasses. In radiometric terms it is the normalized 
bidirectional reflectivity. The diffuse reflectivity (Bond albedo) can now 
be expressed in terms of the normal albedo, namely 

where the Integral is referred to as the phase integral. 

Retroref lectivlty (or geometric albedo ) is the value of the bi- 
directional reflectivity when the influx and efflux angles are equal (zero 
phase angle), but not necessarily zero themselves, namely 

Note that the normal reflectivity is a special case of the retroref lectivity . 
The geometric albedo accounts for the variation in apparent brightness 
across a planetary disk near opposition. Because the lunar disk appears to be 
unlforui, l.e. no limb darkening, its retroref lectivlty is constant and 
equal to its normal reflectivity. This is not true for the martian disk 
where limb darkening effects have been observed. Frequently the terms 
normal albedo and geometric albedo are used Interchangeably leading to 
some confusion. The average reflectivity of a planetary disk near 


8 



opposition is an area weighted average of its retroref lectivity or geometric 
albedo but is sometimes identified as the geometric albedo of the planet. 


In many situations the bidirectional reflectivity, Eq. (2), is not 
a function of the four angles but only three where 

g is the phase angle shown in Figure 1.* This allows the bidirectional 
reflectivity to be written as: 


N _ r ( S- 

- - o J 




( 9 ) 


This simplification is a consequence of assumed symmetry properties of the 
material, e.g. that the reflectivity is independent of a rotation about its 
normal . Although most radiometric ("photometric") studies of planetary 
surfaces accept this assumption it may not always be valid--for example if 
wind blown dust were to have a preferential deposition could be 

ambiguous 

' Finally using the more conventional notation, and^^£:', the 
component of radiance due to a planetary surface is given by the expression: 


S t(6,A) (10) 

where S(<^) is the irradlance on the surface at wavelength A and the 


*Note that J fi-- ^ g £ Q, + 
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loss factor due to the transmissivity of the optical path. Note that in 
the presence of a significant atmosphere S will include both solar and sky 
components and that the apparent radiaiKe of the planet will include an 
atmospheric component as well as a surflclal component. For the moon. 

Mercury and Mars (except during dust storms) the surficial component 
dominates the observed radiance . 

Unfortunately one of the other ambiguities in the definition of 

radiometric terms is that factors of fr "come-and-go" . Some authors may have 

a factor of l/ff multiplying the right side of Eq. (10). This factor originates 

because of a difference in the definition of photometric units between the 

English and metric systems. To further complicate the situation, a similar 

difference does not formally exist for radiometric units although many authors 

impose such conditions in order to use operationally measured reflectivities. 

The English units were developed to facilitate the operational measurement of 

reflectivity such that a perfectly reflecting diffuse (Lambertian) surface has 

a bldirection reflectivity of unity. The English unit of luminance (ft-Lambert) 

has been rationalized by a factor of fT in addition to the conversion factor 
2 2 

from m to ft . Consequently the luminous bidirectional reflectivity, Eq. (2), 
is increased by T7' ,vlz: p'= ^ the diffuse reflectivity, Eq. (4), 

becomes 9d =V (i.e. the bidirectional reflectivity of a Lambertian surface is 
numerically equal to its diffuse reflectivity) and the factor of 1 /tt enters 
Eq. (10). Note also that Nr(Gr-»^r^ constant for a Lambertian surface and the 
definition based upon the ratio of reflected to incident power, Eq. (1) also 
becomes P = ITN^/H^ and the two definitions of luminous reflectivity are identi- 
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cal for Lambertian surfaces. Operationally we frequently measure the bidirectional 
reflectivity of a material by comparing the received power to that received from 
a Lambertian standard (withpd -1.0) at the same orientation, the value obtained 
is governed by the alternate definition discussed here and the multiplying 
factors discussed above apply. 
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SECTION 3 


BIDIRECTIONAL REFLECTIVITY (PHOTOMETRIC) FUNCTION 

It is clear from the results of the previous section that some 
knowledge of the bidirectional reflectivity characteristics of the surficial 
materials on a planet are required in order to interpret the radiance 
measured by either an imaging system or a radiometer. Even if a planetary 
surface were composed of a single type of material with a well defined 
particle size regime, the radiance of the planetary surface will vary due 
to geometrical considerations alone. Both the normal reflectivity and the 
normalized bidirectional reflectivity (photometric) function can change due 
to compositional differences on the planetary surface. Consequently, it 
becomes a matter of removing the geometrical dependence of the radiance 
to arrive at an "albedo" which contains information about the composition 
of the surface. One approach is to use Eq . (10) to convert the observed 
radiance into the normal albedo. This requires knowledge of the normalized 
bidirectional reflectivity or photometric function. The normal albedo 
differences subsequently derived can be due to chemical or mineralogical 
differences, particle size differences or a combination of these effects. 

A primary objective of this study was to determine whether observed albedo 
differences (dark and light markings) in Mariner VI and VII photography of the 
Martian surface were due to chemical effects or particle size differences. 
However the analyses methodology did not require the computation of the 
normal albedo or reflectivity. 
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If we assume that there are no chemical or composition differences 
then the photometric function or normalized bidirectional reflectivity can 
be derived from measurements of the radiance of the surface at different 
Illumination and viewing conditions. The measured data is fit to an analytic 
model for the photometric function. The most frequently used functional 
representation for surface radiance Is the Minnaert equation (Mlnnaert, 1941). 
This equation Is a parametric expression developed to obey the Helmholtz 
reciprocity law -- the bidirectional reflectivity is invariant upon reversability 
of the incident and emission angles, namely 

The resulting parametric equation has the form 




) (C 06 C) 


( 12 ) 


where the exponential parameter, k. Is a function of both the phase angle 
and wavelength and can be thought of as a limb darkening parameter. The 
moon, which has no limb darkening effect, has a value of k=0.5 at zero phase 
while a Lambertian surface would have a value of k=1.0. The Mlnnaert law is 
frequently written with the symbol "B” in place of the symbol "N" which we 
have employed. The former Implies the brightness of the surface in contrast 
to the radiance and consequently In keeping with the discussion in Section 3 
we have avoided its use here. The term Np(g,]^), frequently written as "Bq" 
can be a function of phase angle without violating the reciprocity requirement 
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since It does not depend on the incident or emission angle. At zero phase 
note that Np is equal to the radiance obtained at normal Incidence and normal 
emission, namely (0,^)=N^(0,0)=HQ(i\)|' where is the irradiance at normal 
Incidence. This is required so that Eq . (12) has the correct limit at 
L = 0. The Minnaert equation Is frequently used in a modified form 

obtained by implicitly dividing both sides by Hq(^) . In this case the 
value of the parameter Np or Bq at zero phase angle is the normal 
reflectivity or albedo. -If we assume that Np is not a function of the phase 
angle, then Np:=Hp(;^)p for the reasons given above. Some authors make this 
assumption without noting that it may not be justified. 


Comparing Eq . (12) to Eq . (10) and assuming that the loss due to 
the atmosphere is negligible, the photometric function described by the 
Minnaert equation is 








Hcl,X) C (A) 


V. ) (^0 S (. ^ 


(13a) 


Equating the photometric function and the Minnaert equation, the practice 
of some authors, is somewhat imprecise and misleading. If Np is not a 
function of the phase angle the photometric function becomes 






(13b) 
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The Far Encounter photography obtained during the Mariner VI and VII missions 

has been employed to estimate the parameters Np,k in Equations (12) and (13) 
(Young, 1971). The video signal was converted to reflectivity using available 
calibration data such that the value of Np/Ho at zero phase would be the normal 
reflectivity (albedo), If It could have been estimated. Log-log plots of 
N cose versus cos i- cost called •'Mlnnaert plots" were constructed to 
estimate Np/H^ and k. Since these photographs were obtained with the long 
focal length TV camera and occurred at limited phase angles , 23° "tz® this 

evaluation did not Include the variation of Np/H^ or k with phase angle or 
wavelength. The results obtained from each of the missions are summarized 
In the table below. The variation of k from place to place on Mars was 
attributed to a variation In composition and supported the hypothesis that the 
Martian surface Is more than a two-component system. 


i 

Table 1 

Mlnnaert Parameters from Mariner VI and VII Imagery 


Region 

West 

Longitude, 

deg. 

Latitude, 

deg. 

^6 



N A 
H 

o 

Ster ^ 

N 7 

_eZ 

H 

o 

Ster ^ 

Ophlr 

68 

-8 


m 



Center of Elysium 

213 

+12 

0.56 




Aeolls 

213 

0 

0.61 

mm 



Center of Syrtls Major 

290 

+5 

0.46 

0.48 


0.071 

Soils Lacus 

90 

-40 

0.66 

0.60 

0.117 

0.093 
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Earth-based observations of the radiance of Mars have also been 


employed to estimate the parameter k (Binder, 1972). iUnf ortunately 
Earth-based observations are constrained In that measurements cannot be 
made of the Martian surface at phase angles greater than 40? due to 
ephemeris considerations. The results of these measurements show a slightly 
higher value of k at a 23° phase angle, namely k-^0.8. These data were 
measured at a wavelength of 0.60 ^m. The differences between these two 
measurements of k are most likely the result of the differences In the 
spectral bandpass of the sensors. This study also evaluated the variation 
of k with wavelength with the resulting dependence at approximately 10° 
phase angle . 

k( 0.91 W 

since the spectral passband of the B camera varies from 0.48^ to approximately 
0.65^m (Danielson, 1971) the value of the k parameter according to Eq. (14) 
corrected for the phase angle difference would vary from 0.74 to 0.84, still 
somewhat higher than the Mariner derived values. We note, however, that 
Binder arbitrarily assumed that k at zero phase angle was 0.5 (no limb 
darkening effect) . 

A more recent evaluation of the Mlnnaert coefficients derived from 

observations on the Mariner IX Imagery (Thorpe, 1973) resulted In values 

presented in Table 2 as a function of phase angle at an effective wavelength 

of 0.56 ^(m. In this case also the value of Np at zero phase would be the 

normal reflectivity. Note that the results presented In the table clearly 

Indicate that Np Is dependent on the phase angle. If we fit results of 
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Thorpe for the variation of k with phase angle with a linear equation we find 
that 

k(g) = 0.63 + 0.0035 g (15) 

since k(0)>0.5 this Indicates a limb darkening effect on Mars . 

Measurements of the Mlnnaert coefficient k at a mean phase angle of 
58° were also made during the Initial portion of the Mariner IX mission while 
a dust storm covered Mars (Masursky et al, 1972). For these measurements 
k increased from 0.93 to 1.12 with Increasing wavelength and Indicate a 
nearly Lambertian surface (k=l .0) . This result would be expected for an 
optically thick atmosphere composed of small particles ("dust") . 

A recent research effort has led to an alternative expression for 
the photometric function (Meador, 1975). Although more complex than the 
Mlnnaert law the Meador equation has the distinct advantage that its 
parameters can be related to the physical characteristics of the reflecting 

I i 

surface, e:g. particle size, single-particle albedo and compactness. The 

i ' 

explicit functional form will not be presented here and the interested 
reader Is referred to the cited reference for details. This expression has 
been compared to the three data sets previously used to derive the parameters 
for the Mlnnaert law as discussed above (Weaver, 1974). Based upon the 
fit of the Meador expression to the experimental observations and laboratory 
data on Colorado Basalt it was concluded that the mean intercenter spacing 
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of adjacent particles is about 4/3 of the mean diameter. The mean diameter 
was concluded to be greater than 225 ;um and not In conflict with other 
indications that the mean diameter of Martian surface particles is about 
400jU.m. 

Regardless of which expression is used to represent the bidirectional 
reflectivity or photometric function in order to analyze the radiance of a 
planetary surface as measured by either an Imaging system or a radiometer, 
the values of the Influx or incident angle, 1, the efflux or emission angle, 
e , and the phase angle, g, must be determined. If the latitude and 
longitude of the sub-solar point are known ) the Incident angle 

at a point having latitude and longitude ,"X ) is given by 

COS i - -f- c-OS<p> Cos 

The sub-solar point is the location on the planetary surface where the line 
from the center of the sun to the center of the planet Intersects the surface. 
Similarly if the sub -spacecraft point is known emission angle 

can be determined by using 

I 

cose - cos ^^^CO£(X-'X (17a) 

As the spacecraft approaches the planet a (parallax) correction may be 
required. The emission angle, given by Eq. (17a) is Increased by ^ . 

■hcLr, S> ~ COS6J (17b) 
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where £. is the value determined from Eq . (17a), Rg is the distance from the 

•I 

spacecraft to the center of the planet and r is the radius of the planet. 

For the Mariner VI and VII Imagery the values of Rg have been computed by 
Davis (1971) and r.-l- 3385km for Mars. For the Mariner VI and VII imagery 
Rg/r~2-3 and the correction factor is significant. The locations of the 
sub-solar and sub-spacecraft paints are not as readily documented and must 
be determined from orbital data. If the azimuth angle between the incident 
and emission planes is known and equal to ^ then the phase angle is given by 

' i c-Ose:- -f Cos& 

^ (18) 

The angle & is the difference between the bearing of the sub-solar point, Zg^, 

Sml co£<p 

and the bearing of the sub-spacecraft point, Zgg, 

Cc5S/ ~ 5/>7^ (19b) 

The sign of Zg 2 and Zgg is determined by the value of the longitudes ^ g 2 ss 
with respect to^ . 
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SECTION 4 


RADIOMETRIC APPLICATIONS IN GEOLOGY 

In this section we discuss some of the applications of radlometry 
and the bidirectional reflectivity function In obtaining information about 
the geology of planets. No attempt has been made to make this discussion 
exhaustive but only representative of some of the applications available. 

The discussion has been divided Into three categories which Include (1) the 
Influence of reflectivity on topographic Information, (2) the use of 
reflectivity for global geology, and (3) the use of local changes 
in reflectivity In assessing con^osltlonal differences in surflclal 
material. 

4.1 Topographic Applications 

Two examples where the reflectivity characteristics of a planetary 
surface influences topographic information include the use of imagery in 
crater distribution studies and the development of topographic profiles from 
photometrically corrected image data. 

The effect of the finite resolution of an imaging sensor (e.g., 
Mariner television systems) upon the accuracy of distribution data is well 
known, namely as the diameter of the crater approaches the resolution limit 
of the imaging system the measured distribution falls below the actual 
distribution due to the inability to identify and measure the craters on the 
reconstructed imagery. The reflectivity characteristics of the planetary 
surface play a role in how rapidly the loss of information due to limited 
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resolution occurs. This point is not necessarily appreciated by all investi- 
gators utilizing the imagery in obtaining their photogrammetric information. 

The resolution of an imaging system is most often quoted assuming that the 
object being imaged is high contrast (having extreme radiance values). As the 
contrast of the object decreases the resolution limit becomes poorer. That is, 
the inherent resolution limit of the imaging system is only achieved under 
Illumination and viewing conditions where the objects under study have a 
reasonably high contrast. At high phase angles most surface features have 
associated shadows which enhance their contrast. Consequently, one would 
expect the Information content of Imagery acquired at large phase angles 
to be higher than that acquired at lower phase angles. As part of the 
development of quality evaluation techniques for Lunar Orbiter photography 
the author (Kinzly, 1967) develope’d an expression for information content 
in terms of a number of mission parameters including phase angle which was 
subsequently related to the minimum detectable crater diameter of Lunar 
Orbiter photography. The variation of the minimum crater diameter with 
respect to phase angle for this photography is shown in Figure 2. A 4;1 
variation occurs as the phase angle changes from 0 to 90°. The change is 
most dramatic in the region from 60° to 90° emphasizing the value of large 
phase angle photography. More recently the importance of illumination 
conditions upon comparative planetary geology has been evaluated (Schultz, 1976). 
Comparative geology uses Imagery from several planetary and lunar missions. 

It is noted, for example, that the information content of the Mariner IX 
Imagery is comparable to the Mariner VI and VII imagery despite the smaller 
slant ranges Involved on the former mission. The improved theoretical 
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Figure 2 EFFECT OF PHASE ANGLE ON INFORMATION CONTENT OF NEAR VERTICAL 
LUNAR ORBITER PHOTOGRAPHY 
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resolution limit of the Mariner IX Imagery Is offset by the smaller phase 
angles at which the Imagery was obtained. A comparison of TV and photographic 
imagery of the moon shows that image enhancement allows the TV Image to 
approach the limiting resolution of 2.2 TV lines. However, the enhancement 
can distort topographic features which are on the order of several resolution 
elements. 


Since the radiance of a point on a planetary surface depends upon 
the incident and! emission angles measured with respect to the local surface 
normal, an accurate measure of surface radiance can be Interpreted In terms 

i 

of the surface orientation if the bidirectional reflectivity function for the 
surface is known. The conversion of surface radiance measured from an Image 

Into local slope and subsequently into a topographic profile has been termed 

} 

"photocllnometr.y" and was applied to both Ranger (Rlndf leisch, 1966) and 

{ . 

Lunar Orblter (Lamblotte, 1967) Imagery. In the Lunar Orbiter application, 
the technique was used to derive relative surface roughness indices at 
potential Apollo landing sites. An evaluation of the quality of Lunar 

Orbiter imagery. (Kinzly, 1968) showed that the medium resolution photographs 

! 

of Missions I and II received excessive exposure and the image densities could 
not be reliably converted to surface radiance. The limitation that this 
placed on the application of photocllnometry was investigated (Cambell, 1968). 
More recently photocllnometry was used on Mariner IV Imagery to produce depth/ 
diameter data for martian craters (Cintala, 1976) . Because of the nonlinear 
relationship between surface radiance and surface slope, the average radiance 
measured by an Imaging system Is not necessarily equivalent to the average 
slope of the surface -- both averaged over a resolution element. Consequently, 
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the resolution of the imaging system can affect the accuracy of the profile 
Information generated by photocllnometry . This analytical procedure, to the 
best of our knowledge, has nevery received widespread application beyond 
Lunar Orblter. 

4 .2 Definition of Global Geology 

As noted In Section 3, differences In the normal albedo of a 
planetary surface can be produced by chemical or mlneralogical differences, 
particle size differences or a combination of these effects. Consequently, 
a geologist frequently makes use of the normal albedo of the surface as an 
aid In defining geologic units. This potential utilization led to the 

development of albedo maps of the moon (Pohn, 1970), Mars (Cutts, 1971 and 

i 

de Vaucouleurs, 1973) and Mercury (Dzurlsln, 1976). In addition to the fact 
that albedo boundaries can be coincident with the boundaries between different 
geological units, correlation between reflectivity of surflclal materials 
has been noted to vary Inversely with age on the lunar surface -- older and 
more heavily cratered areas having a higher reflectivity. 

In addition to the production of albedo maps which are used in 
qualitative geologic evaluations, analysis of the reflectivity characteristics 
of the light and dark areas of a planet on the global scale have been used to 
draw Inferences on the composition of surflclal materials. Typical examples 
of these type of analyses have included those of Pollack and Sagan (1967) 
where It was concluded that the bright and dark areas of Mars appeared to have 
a very similar chemical composition with goethite, a major constituent of both 
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areas. It was also suggested that both areas ijere covered by a fine powder 
and that the seasonal change In the reflectivity of the dark areas was produced 
by aeolian transport of material to and from the dark areas. The average 
diameter of the particles in the bright areas is estimated to be 50 yUm while 

the dark areas were estimated to have particles ranging from 200 ^im to 400 yum 
depending upon the period relative to the seasonal darkening. Based on 
visible polarlmetry and IR radlometry they proposed llmonite as a primary 
constituent of this powder. 


A second example of global reflectivity analysis Is that of Binder 
and Jones (1972) . Observations were made using a ten channel spectroradlometer 
operating from 0.60 - 2.27 They concluded that the measured spectro- 

reflectlvity fell into two well-defined groups, mare and desert units, 
indicating that the martian surface consists of two types of materials. By 
comparison with laboratory data they proposed that both surface units contained 

llthlc soll|s: having a limonlte stain and that the soils of the maria are 

1 ( 

richer in pyroxene, olivine, or both than the desert soils. This investigation, 

1 I- 

therefore, suggests that the differences between martian maria and deserts 
are due to composition and not due to differences in particle size as proposed 
by Pollack and Sagan. 

A third example of global reflectivity investigations is the work 
of McCord and his associates (McCord and Adams, 1968; Adams and McCord, 1969; 
McCord and Westphal, 1971). These investigations involved spectroradiometric 
observation of Mars from 0.30 - 2.5 ^m. Comparison of the resulting spectro- 
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reflectivity curves to laboratory data led to the conclusion that oxidized 
basalts are a significant constituent of the martian surface. Evidence has 
been suggested for the existence of basaltic flows and mare on Mars, 
thereby supporting this hypothesis. It should be noted, however, that the 
Mariner IX Imagery acquired after a martian dust storm contained numerous 
bright and dark areas produced by dynamic aeollan transport of material over 
the martian surface. These observations support the hypothesis by Pollack 
and Sagan that the seasonal changes In reflectivity are a result of a change 
in the average particle size. It Is, of course, quite possible that both 
mechanisms play a role in the dynamic characteristics of the reflectivity of 
the martian surface. 

I 

4 .3 Interpretation of Local Reflectivity Differences 

The evaluation of local changes in planetary surface reflectivity 
requires the availability of radloraetrically calibrated high resolution 
imagery. Consequently, the techniques available for such analyses are less 
developed although they parallel those used In evaluation of reflectivity 
changes on a global scale. The distinction that we make between global and 
localized analyses can be correlated directly with earth-based versus space- 
craft observations of the planet. An example of the qualitative analyses of 
local reflectivity variations is the use of the radial bright ray patterns 
associated with some craters on the lunar surface in order to derive relative 
ages. Quantitative techniques are available as illustrated In this report 
to analyze the dark and light markings occurring In the Mariner VI and VII 
Imagery . 
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As noted In the previous section. Young and Collins (1971) analyzed 
the Mariner VI and VII far encounter photographs to derive the parameters of 
the Mlnnaert photometric function for five regions of the martian surface. 

We consider this analysis to be global rather than local, however, since it 
only employed the far encounter pictures wiilch have resolutions ranging from 
10 - 100 km. 


A method for evaluation of local changes in the reflectivity of 
soils employing radiometric analysis using conventional color aerial photography 
has been recently reported (Plech, 1974). The technique was developed to 
assist in terrestrial soil surveys by supplementing conventional land form 
analyses with reflectivity information extracted from color imagery. The 
cause of reflectivity variations in a soil unit are evaluated by computing the 
ratio of the reflectance in two spectral bands. Information on the relative 
soil moisture and texture characteristics is derived from the reflectivity 
ratios obtained from the imagery. In this case the aerial camera is employed 
as a radiometer. In addition to calibration^ of the object -to-lmage radiometric 
response^ the effects of the atmosphere must be removed from the imagery in 
the earth-based application. This is accomplished by radiometric analysis 
of standard scene objects such as shadow areas. If the darker soil element 
has a greater red-to-blue reflectivity ratio than the lighter soil element, 
the reflectivity variation is caused by differences in moisture content. If 
the darker soil element has a smaller ratio, the decrease in reflectivity is 
produced by an increase in the average particle size. Rather than computing 
the reflectivity ratios in two spectral bands the equivalent procedure of 
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computing the ratio of the reflectivity of the darker unit to the lighter unit 
in each of the spectral bands selected can be used. If the ratio of the 
darker to lighter soil element at the longer wavelength (e.g. red) is less 
than the ratio at a shorter wavelength (e.g. blue) then the change in 
reflectivity can be attributed to a larger particle size for the darker 
element. Figure 3 shows laboratory data obtained in support of this analysis 
technique. The figure shows the ratio of the darker to lighter (i.e. wet to 
dry) soil elements for several’ terrestrial sands and varying amounts of 
moisture content. The reflectivity ratios were obtained from densitometer 
measurements of vertical photography (normal emission) obtained under ambient 
sunlight such that the incident angle and phase angle varied from 30° to 50°. 
Note that ^he ratio Increases with wavelength supporting the interpretation 
of reflectivity ratios proposed above. In the case of reflectivity changes 
due to particle size variation, the ratio decreases with Increasing wavelength. 
Additional experimental data is presented in the cited reference. This size- 
reflectivity change rule is not applicable to all types of minerals that 
could be constituents of soils. The reflectivity versus size effect has 
been investigated for likely martian surface materials (Salisbury, 1968) 
and is discussed further in Section 5.2. 

i I 

I 

The application of this analysis technique to adjacent light and 
dark surflclal elements on the martian surface can be used to distinguish 
between a chemical cause for lower reflectivity, namely absorbed moisture or 
a physical cause due to varying particle size. Note that if the composition 
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Figures WET/DRY RATIOS DETERMINED PHOTOGRAPHICALLY FOR VARIOUS 
TERRESTRIAL SOILS AND MOISTURE CONTENTS BY WEIGHT 
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of the light and dark areas Is assumed to be similar or if the bidirectional 
reflectivity is assumed constant, the ratio of the radiance of the dark to 
the light area is equal to the ratio of their reflectivities and no explicit 
correction for the bidirectional reflectivity or photometric function is 
required. This technique was applied to the analysis of adjacent light and 
dark areas occurring in the Mariner VI and VII imagery -- the results are 
presented in the next section of this report. 
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SECTION 5 


ANALYSIS OF MARINER VI AND VII IMAGERY 

This section of our report documents analyses carried out on 
Mariner VI and VII Imagery to Interpret local variations In reflectivity by 
using reflectivity ratios derived from the Imagery. Darkening of the martian 
surflclal material can be produced by chemical effects such as absorbed water, 
by changes In the particle size or by changes In the composition of the 
surface material. Identification of the presence of significant amounts of 
moisture In the martian soil Is important to the assessment of life forms on 
Mars while the assessment of relative particle size differences Is Important 
to the evaluation of dynamic aeollan processes. Alternate near encounter 
frames from the wide angle Mariner cameras were taken through different 
spectral filters (red, green or blue) and the overlap between successive 
frames show the same area of the martian surface in two spectral bands » These 

I 

regions of [overlap provide an opportunity to use reflectivity ratio analysis 
techniques In the study of martian planetology. In this section we (1) 
Identify the areas where bl-band spectral data is available, (2) present 
additional data in support of the interpretation of the reflectivity ratios 
beyond that identified previously in Section 4.3, (3) describe the tools 
developed to compute the reflectivity ratios from digital data supplied 
by the Image Processing Laboratory at JPL and (4) present the results of 
the analysis of adjacent light and dark areas occurring in the Mariner 
photography . 
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5.1 Mariner VI and VII Multi-Spectral Image Data 

The objective of this Investigation was to Illustrate the utility 
of radiometric analysis, particularly reflectivity ratios, to planetology in 
general and the study of Mars in particular. We planned to use TV Image data 
obtained from the Mariner VI, VII and IX missions. Mariner VI and VII were 
fly-by missions that had a variety of experiments Including a dual TV Imaging 
system (Danielson, 1971) . One TV sensor employed a short focal length lens 
(■^SOmm) while the second sensor employed a long focal length lens (^500mm) 
and were termed the "wide angle" and "narrow angle" cameras respectively. 

Since the imagery was obtained during fly-bys the resolution (^2.2 TV lines) 
changed from image to image. For the imagery utilized It was approximately 
2-4 km for Mariner VI and 4-5 km for Mariner VII. The wide angle camera 
contained a filter wheel with four spectral filters. The spectral transmission 
of the filters incorporated into Mariner VI and VII are shown in Figure 4. 
Alternate near encounter frames from the wide angle Mariner cameras overlapped 
such that the same area of the martian surface was taken through two 
different spectral filters. The Initial step in the research effort was the 
identification of those areas where bi-band coverage was available. 

Reflectivity ratios computed within these areas could be used for geologic 
mapping, discrimination of atmospheric phenomena, obtaining information about 
polar cap structure and obtaining information about the cause of local 
variations In the reflectivity of the martian surficial material. 

Six areas of bi-band coverage were identified from the Mariner VI 
mission and six areas of bl-band coverage from the Mariner VII mission. 
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Figure 4 SPECTRAL TRANSMISSION OF MARINER VI AND VII WIDE ANGLE CAMERAS 
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Table 3 summarizes these areas which are designated by the mission number 
followed by a letter code A through G. This table Identifies the location 
of each of these areas and the frame numbers and spectral filters of each 
of the overlapping photographs. Figure 5 shows the location of each of the 
bi-band coverage areas on the 1971 Shaded Relief Map of Mars. 

A number of potential bi-band coverage areas were eliminated from 
consideration. These Included the two areas of overlap between 7N5, 7N7 and 
7N9 . These photographs were obtained at a highly oblique angle and 
had lower resolution thau the subsequent near encounter frames. Consequently 
the radiometrlcally corrected digital imagery was not requested from JPL. 

In retrospect this is somewhat unfortunate since these photographs included 
some of the areas in Meridian! Sinus coincident with areas 6B and 6C previously 
selected. In addition, the overlap area between frames 7N29 and 7N31 

were deleted from consideration because of their proximity to the terminator 

I ’ 

which yields a lower quality due to the low exposure level on the vidicpn. 

One area of overlapping coverage among frames 7N15, 7N17 and 7N19 
actually had tri-band coverage. This is a small area in the South Polar Cap 
located at 72°S - 80°S, 332°W - 348^^7. Although no specific analysis was 
undertaken over this area, several features appear that might be caused by 
atmospheric phemonena. Figures 6 through 9 show the Mariner VI bi-band 
coverage areas on enhanced versions of the imagery transformed to othographic 
projections. Figures 10 through 13 show the bi-band coverage areas of 
Mariner VII on similar projections. Digital tapes containing TV image data 
for each frame identified in Table 3 were obtained from Mr. A. Collins at the 
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Table 3a 

Mariner Bl-band Coverage Areas 
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Mariner Bl-band Coverage Areas 
Mariner VII 
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Cratered Deposits, Undivided 
Sparsely Cratered Plains 
Mountainous Deposits 






W. LONGITUDE 


Figure 6 BI-BAND AREA 6A SHOWN ON ORTHOGRAPHIC PROJECTION OF FRAME 6N7 


. LATITUDE S. LATITUDE 





W. LONGITUDE 



Figure 7 BI-BAND AREAS 6B AND 6C SHOWN ON ORTHOGRAPHIC PROJECTION OF FRAME 6N11 


40 


L LATITUDE S. LATITUDE 


W. LONGITUDE 



Figure 8 BI-BAND AREAS 6D AND 6E SHOWN ON ORTHOGRAPHIC PROJECTION OF FRAME 6N17 
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S. LATITUDE 









W. LONGITUDE 



Figure 10 BI-BAND AREAS 7A AND 7B SHOWN ON ORTHOGRAPHIC PROJECTION OF 
FRAME 7N13 
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LATITUDE 


W. LONGITUDE 

360 ° 350 ° 340 ° 330 ° 320 ° 310 ° 



Figure 11 BI-BAND AREAS 7B, 7C AND 7D SHOWN ON ORTHOGRAPHIC PROJECTION OF 
FRAME 7N17 
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S. LATITUDE 



W. LONGITUDE 

336 ° 332 ° 328 ° 324 ° 320 ° . 316 ° 312 ° 



Figure 12 BI-BAIMD AREAS 7E AND 7F SHOWN ON ORTHOGRAPHIC PROJECTION OF 
FRAME 7N25 
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LATITUDE 



W. LONGITUDE 



Figure 13 BI-BAND AREAS 7F AND 7G SHOWN ON ORTHOGRAPHIC PROJECTION OF 
FRAME 7N27 
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S. LATITUDE 




Jet Propulsion Laboratory. The data obtained was the LMICOR version which 
was produced by the Image Processing Laboratory at JPL by correcting the raw 
TV data for radiometric distortion introduced by the TV systems. The corrected 
data represents the surface bidirectional reflectivity multiplied by a factor 
of 6.0 in order to optimally utilize an 8-blt foirmat. 

The area that received the most detailed analysis was 6C . This 
area has light and dark markings on crater floors which were evaluated to 
determine whether this variation is caused by absorbed water or by changes 
in particle size. The method employed uses the reflectivity ratios in the 
manner discussed in Section 4.3, Section 5.2 contains additional data on the 
relationship between particle size and reflectivity that supports our analysis 
of these markings. The results of the analysis Itself are contained in 
Section 5.4. 

5.2 Reflectivity - Particle Size Effect 

A size-reflectivity change rule was described in Section 4.3 that 
has been utilized in terrestrial soil surveys. This rule indicated that the 
reflectivity ratio of the dark-to-llght soil element should decrease with 
increasing wavelength if a change in particle size was responsible for the 
change in reflectivity. 

The proposal that llmonlte be a major constituent of the martian 
surface led Salisbury and Hunt (1968) to conduct a study of the spectral 
behavior of likely martian surface materials. They classified potential 
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materials by their change in reflectivity as a function of particle size. 

The class of transparent materials, which includes silicate minerals, show 
a Increase in reflectivity with a decrease in particle size as 
assumed in the development of the size-reflectivity change rule referenced 
above. Opaque material, particularly metal sulfides, showing a decrease in 
reflectivity with decreasing particle size. A third class of materials which 
Salisbury and Hunt term "trans -opaque" exhibits a decrease in reflectivity 
in the opaque portion of the visible spectrum and an increase in reflectivity 
in the transparent portion of their visible spectrum. That is, the spectral 
reflectivity curve for one particle size will cross that for a different 
particle size leading to a reversal in reflectivity between the two particle 
sizes within the visible spectrum. This behavior is exhibited by several 
ferric oxides such as limonlte, goethite and hematite that are likely 
candidates for existence on the martian surface. As Salisbury and Hunt point 
out, light and dark markings on the martian surface composed of high concen- 
trations of these types of materials should exhibit a contrast reversal at 
some point in the visual spectral region. That is, the light markings should 
become dark and the dark markings become light as the radiometric or image 
data Increase in wavelength. Since such a reversal has not been observed 
Salisbury and Hunt concluded that the hypothesis that the martian soil 
consists in a large part of limonlte and that the reflectivity differences 
are due to particle size changes are Incompatible. As an alternative they 
propose that the most likely soil is one composed of silicates lightly stained 
or coated with ferric oxides . In terms of the present evaluation of the 
Mariner Imagery the dark-to-light reflectivity ratio could Increase with 
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wavelength between the blue and green images and subsequently decrease with 
wavelength between the green and red Images If a significant amount of the 
trans-opaque type materials composed the elements under evaluation and the 
change in particle size was the principal cause for the change in reflectivity. 
Our analysis concentrated on the local reflectivity changes occurring in 
bl-band coverage area 6C involving overlapping green and red Imagery. 

To further support our analyses using the reflectivity ratio, the 
ratios for typical martian materials were evaluated over the spectral responses 
of the Mariner VI and Mariner VII wide angle cameras presented earlier in 
Figure 4. Hunt and Salisbury (1970, 1971) have obtained a number of spectral 
reflectivity curves for various rocks and minerals as part of a laboratory 
study of spectroscopic remote sensing techniques. We selected typical candi- 
dates from the published data for silicate minerals (1970) and oxides and 
hydroxides (1971) . The spectral reflectivity for each of the minerals selected 
is shown in Figure 14 (although measurements were made at IR wavelengths this 
data was not of interest here) . The reflectivity data were obtained near 
normal incidence and emission angles at a phase angle of 15° and referenced to 
a Lambertian surface with a diffuse reflectivity near 1.0. Generally the data 
were obtained for four different particle size regions allowing six different 
ratios of the reflectivity of larger-to-smaller particle sizes to be obtained. 
For materials whose reflectivity increases with decreasing particle size this 
ratio is equivalent to that of the darker- to-lighter soil elements. The six 
different ratios obtained are identified in Table 4 using the letter s}nnbol A 
through F. 
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Figure 14 SPECTRAL REFLECTIVITY FOR TYPICAL MARTI AN MINERALS AS A 
FUNCTION OF PARTICLE SIZE 


50 








Table 4 


Large-to-Small Particle Index 


1 

Symbol 

Particle Size Ratio 

A 

250-1200 Am/0-5>tm 

B 

250-1200/^('m/0-74xm 

C 

250-1200//m/74-250/<'m 

D 

74-250^4( m/0-5 

E 

74-250"a a/0-74/<m 

F 

0-74 /< m/0-5 ^m 

^ 


From the silicate minerals olivine and p5rroxene were selected based 
upon the analysis of Binder and Jones (1972). Montmorlllonite was arbitrarily 
added as the third selection from the silicate minerals. Limonite was 
selected from the ferric oxides since It has been proposed as a major 
constituent by Pollack and Sagan (1967) and as a stain on silicates by 
Binder and Jones (1972) and Salisbury and Hunt (1968). Pollack and Sagan 
proposed goethite as a major constituent of the martian surface and con- 
sequently this was also selected by the oxide group. Finally, Hematite was 
added as a third ferric oxlde--a logical addition to those already selected. 
Table 5 presents the average reflectivity of these minerals within the 
Mariner VI and VII spectral bands. As the table shows, the reflectivity 
increases with decreasing particle size In all cases except for the smallest 
particle size group of goethite In the green band (and presumably the blue 
band) and for all of the particle size groups measured In the case of 
limonite for both the blue and green spectral bands. Consequently, the 
ratio of the larger -to-smaller particle size groups obtained from this data 
are equivalent to the ratio of the darker -to-llghter surface elements except 


for the cases cited. 


51 






Table 5 


Average Reflectivity of Potential 
Martian Materials Within the Mariner Spectral Bands 


Mineral 


0-8 Goethlte 


Size Range 


- 1200 u 

- 250 

- 74 u 

- 5 // 


Mariner VI 


Blue Green 


Mariner VII 


0-9 Hematite 

74- 250 
0- 74 ,^m 

0- 5 (<m 

0-11 Llmonlte 

250-1200 yfm 
0- 74//m 

0 - 5 m 

S-llB Mont- 
morlllonlte 

250-1200 //m 
74- 250/^m 
0- 74/(m 

0- 5/|m 

S-14B Olivine 

250-1200// m 
74- 250/4 m 
0- 74/4m 

j 0- 5//m 

S-17C 

Pyroxene 

1 250-1200y4m 
! 74- 250// m 
' 0- 74//m 

0- 5/<m 
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Figure 15 shows the reflectivity ratios determined for the ferric 
oxide samples and the Mariner VI spectral response functions. The ratios were 
plotted at the effective focal length for each of the spectral filters as 
specified by Danielson (1971) . Figure 16 is the corresponding plot for the 
silicate samples. The results obtained In the case of Mariner VII are 
similar to the reflectivity ratios shown here. We would expect that the 
ratios would approach unity as the particle size differences between the two 
groups being ratloed decreases. In addition, the rate of change of the ratio with 
respect to wavelength should decrease as the particle size differences decrease, 
i.e., the slope of the curve becomes less. If we compare the progression of 
curves for cases A to B to C we would expect that these phenomena would be 
apparent In Figures 15 and 16. Except for the cases where a reversal In the 
change of reflectivity with particle size occurs, this phenomena Is exhibited 
for all of the minerals evaluated. A similar trend Is exhibited in going from 
the curves with index D to the curves of index E. We can infer from these 
results that information about the relative dlffersnce in particle sizes can 
be obtained by examining the level and slope of the reflectivity ratio versus 
wavelength curve. The difference between the ratio from the green to the red 
band is presented in Table 6. In this table a positive number indicates a 
decrease in the ratio as the wavelength is increased as expected from the 
size-reflectivity change rule stated previously. Note that most of the 
changes are positive (i.e., show a decrease with increasing wavelength) 
except for olivine which shows an increase in the ratio for all cases. If 
we eliminate this sample from consideration the average value of the change 
which is also presented in the table is greatest for the A ratio which has 
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LARGER-TO-SMALLER PARTICLE SIZE REFLECTIVITY RATIO 





WAVELENGTH (Mm) 


Figure 15 REFLECTIVITY RATIOS FOR SELECTED OXIDES AND VARIOUS PARTICLE 
SIZE GROUPS 
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Figure 16 REFLECTIVITY RATIOS FOR SELECTED SILICATES AND VARIOUS PARTICLE 
SIZE GROUPS 
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Green-to-Red Decrease in the Reflectivity Ratio 
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the Largest particle size differences and least (negative) for the C ratio 
where the particle size differences are not as great. Consequently, we 
conclude that the slope of the reflectivity ratio versus wavelength curve 
Indeed does contain Information about the relative particle size of the 
materials producing the light and dark markings. 

5.3 Evaluation of the Light and Dark Markings Associated with Craters 

Bl-band coverage area 6C contained several craters with light and 
dark markings located Interior to the crater on their floor. Several of these 
craters were chosen for subsequent analysis of the reflectivity ratio of these 
markings from frames 6N11 and 6N13 . The craters selected are shown In 
Figure 17. The light and dark markings were located by sample and line 
number and subsequently retrieved from the LMICOR Images contained on digital 
tape supplied by JPL . Software programs were developed to retrieve this data 
over a rectangular area s'^-eclfled by the Initial and final sample and line 
number. A description of this software Is presented In the appendix to this 
report. Table 7 contains the average reflectivity and subsequent ratio for 
each of the selected dark and light markings In the red and green spectral 
bands . 


These data are not the only source of reflectivity ratios of dark - 
to-llght markings on the martlna surface. Dollfus measured the dark--to-llght 
area contrasts In several wavelength bands. These data were employed by 
Pollack and Sagan (1967) and were extracted for comparison to our results. 
These results are presented In Table 8. Other studies utilizing reflectivity 
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Figure 17 CRATERS ANALYZED IN BI-BAND COVERAGE AREA 6C 
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LATITUDE 


Table 7 


Reflectivity and Reflectivity Ratio of Crater Markings 


Target 

' ' ■ ■ ■ ■ — . - 

6N11 -Green ( .53yVm) 

6N13- 

■Red (.58^m) 

Dark 

Light 

Ratio 

Dark 

Light 

Ratio 




.977 

.0742 

.0879 

.844 

2 - Crater 

KSS 


.921 

.0708 

.0843 

.840 

3 - Crater 

.0592 


.925 

.0658 

.0816 

.806 

4 - Crater 

.0604 

.0632 

.956 

.0668 


.886 

5 - Crater 

.0598 

.0641 

.933 

.0665 


.867 

6 - Crater 

.0603 

.0651 

.926 

.0717 

.0806 

.890 


.0599 

.0638 

.939 

.0734 

.0764 

.961 

Average 

.0595 

.0634 

.940 

.0699 

.0804 

.871 


.0009 

.0015 

.020 

,, ■■ — ■ 

.0035 

.0046 

.049 


A= Green-to-Red Change = .07 


Table 8 

Normal Reflectivity and Dark-to-Llght Reflectivity 
Ratios Derived from Dollfus Contrast Data 



Dark 

Areas 

Bright 

Areas 

Ratio 

.45 

.065 

.071 

.92 

.50 

.098 

.120 

.82 

.55 

.120 

.164 

.73 

.60 

.139 

.210 

.66 

.65 

.150 

.250 

.60 


= Effective Green-to-Red Change = .09 
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ratios were conducted by McCord (1969) and by Cutts (1971). McCord studied 
selected dark and bright regions primarily the Arabia-Syrtus Major pair at 
21 narrow spectral passbands In the visible spectrum. From his results, 

McCord concluded that the bright regions are much redder than the dark 
regions, that Is the ratio of the reflectivity of the dark-to-llght regions 
in the red region Is greater than the blue region. According to the size- 
reflectivity change rule, this would indicate that particle size effects are 
probably responsible for the differences between the bright and dark martian 
regions. Cutts subsequently employed the McCord data to augment his evaluation 
which was obtained using the late far encounter images obtained during the 
Mariner VII mission. Cutts evaluated nine different regional areas ranging 
from high to low reflectivity. The reflectivity ratio for dark- to-light areas 
for selected data presented by Cutts is contained in Table 9. The first 
three ratios were selected to represent that of the darkest- to- lightest 
features. The second two ratios were selected because they represent moderately 
dark-to-llght features. Note that the value of the reflectivity ratio is 
closer to one for the second group than the first group and that the change 
or average slope of the reflectivity ratio curve is less for the second group 
compared to the first group. The third group contains two cases, the first 
being the ratio between two relatively dark features and the second being the 
ratio between two relatively light features. Note in this case that the 
reflectivity ratios are much closer to unity and that the average slope is 
less than either of the two proceeding groups. 
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Table 9 


Dark-to-Light Reflectivity Ratios Derived from 
Selected Cutts/McCord Data Within Mariner Spectral Bands 



Description 

Blue 
( .47/< m) 

Green 
( .53/<m) 

Red 

( .58^ m) 

(Green-to-Red) 

H 

Syrtis Major /Arabia* 

.87 

n 

.65 

.12 

■ 

Margarltifer Slnus/Moab 

.83 

Bn 

.59 

.14 

1 

Meridian! SinusA'Ioab 

.84 

mm 

.65 

.10 

II 

Deucallonls Reglo/Moab 

.97 

.95 

.93 

.02 


Thymlamata/Moab 

.92 

.87 

.80 

.07 

III 

Meridlanl Slnus/Sabaeus 

.95 

.94 

.90 

.04 


Sinus 





1 

Edotn/Moab 

1.01 

1.00 

.97 

.03 


*After Cutts 


All three data sources were plotted for comparison to the Mariner VI 
crater markings in Figures 18 and 19. Figure 18 shows that the crater markings 
show a change in ratio which indicates that particle size differences are 
responsible for the change in reflectivity (the ratio decreases with increasing 
wavelength). Furthermore, comparison of the Mariner results to that of Cutts 
and McCord indicates that the Mariner ratios represent more moderate changes 
in particle size than that due to the extreme reflectivity changes. The 
change in reflectivity ratio from green-to-red of 0.07 is compatible with those 
in Group II in Table 9 and is consistent with the changes obtained in Section 4.2 
for particle size groups B and E. A similar result is exhibited in Figure 19 
when compared with the Dollfus ratios. The Dollfus results appear to represent 
a larger change in reflectivity than that due to the crater markings. The 
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Figure 18 COMPARISON OF MARINER REFLECTIVITY RATIOS TO THOSE OF 
CUTTS AND McCORD 
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Group I data of Cutts and McCord (Table 9) and that of Doll f us (Table 8) are 
consistent with size Groups A and D In Section 4.2. Note that the size ranges 
represented by Groups B and E are compatible with the postulated mean particle 
sizes for the light and dark areas offered by Pollack and Sagan (1967). 

If we compare the reflectivity presented In Table 5 to an extrapolation 
of the Mlnnaert parameter Ng obtained by Thorpe (Table 2) we find that the 
normal reflectivity of olivine appears to be a factor of 2-3 greater than the 
average normal albedo. Consequently, llmonlte stained olivine Is not considered 
to be a likely candidate as a major constituent on this local area. 

Llmonlte stained pyroxene appears to be the best candidate based upon the data 
obtained during the course of this effort. Of course this conclusion must be 
somewhat tenuous since a number of mineral compositions could yield the results 
obtained . 

5 .4 Mariner IX Imagery and Other Related Data 

We had hoped to continue the present effort using Mariner IX Image 
data obtained through multiple spectral filters. It was anticipated that a 
large number of multl-spectral coverages would be available. Unfortunately, 
the filter wheel mechanism became Inoperative during P.evolutlon 118 and only 
a very limited amount of multl-spectral data was obtained. For the subsequent 
portion of the Mariner mission Image data was obtained using the orange 
polarizing filter contained in the fifth filter wheel position. During the 
Initial part of the mission the martian surface was obscured by a dust storm 
and consequently no useful image data of surface features were obtained in 
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this part of the mission when the filter wheel was operative. Mariner IX 
did confirm that dynamic aeollan processes play a significant role In the 
reflectivity variations on the martian surface and that a significant number 
of these variations are caused by changes In particle size. 

An analysis of the light and dark markings present In the Mariner IX 
Imagery obtained subsequent to the subsiding of the dust storm was made by 
Arvidson (1974). As a result of his analysis he found that "dark splotched 
craters In regions with bright streaks usually have upwind bright patches, 
suggesting that these features formed by dumping of bright dust over crater 
rims ... This characteristic pattern minus the bright streaks downwind 
of the craters are identical to those that we have analyzed. Furthermore, 
they have a North to South orientation consistent with the direction of the 
wind pattern observed subsequent to the dust storm through analysis of the 
wind-blown streaks. In addition, a sequence of wind tunnel experiments have 
been carried out by Greeley and his associates (1974). The results of these 
experiments have generated light and dark markings similar to those observed 
during our effort and have further demonstrated that some of the dark markings 
occur from wind erosion and that some of the light markings are deposltlonal 
In nature. The reader la also referred to a comparison of Mariner VI and VII 
imagery and Mariner IX Imagery which Includes the bl-band coverage area we 
analyzed for a qualitative assessment of effect of the dust storm on local 
reflectivity variations (Ververka, 1974). 
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Though it was Intended to demonstrate the use of reflectivity ratios 
as an aid to geologic mapping, the lack of availability of Mariner IX imagery 
greatly degraded this objective. The Mariner VI and VII imagery within the 
bi-band coverage areas occurred over a limited number of geologic units and 
the detailed geologic maps for the Quadrangles containing these areas had 
not as yet been produced. It was also Intended to obtain information about 
the Kiartlan atmosphere through radiometric measurements of the shadows 
associated with geological features, especially bowl-shaped craters. Although 
a number of shadows were Identified in the narrow field of view imagery from 
Mariner VI, in particular Frame 6N22, the lack of available shadows is a 
direct result of the large diameter /depth ratios observed for Mars compared 
to Mercury and the Moon (Clntala, 1976). The diameter to depth ratio for 
craters near 10 km in diameter is approximately 8 times larger while for 
craters near 100 km in diameter it is approximately 3 times larger . The lack 
of Mariner IX bl-band image data further inhibited the evaluation of 
atmospheric phenomena. 
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APPENDIX 


ANALYSIS SOFTWARE 

This appendix describes three items of analysis software which were 
prepared during the course of this study. The principal software package 
consisted of the main program written in Fortran language and an associated 
subroutine, EPIC, written in Assembly language. Listings of the main program 
and subroutines programs are included in Figures A-1 and A-2 respectively. The 
purpose of this program is to retrieve a rectangular array of reflectance 
data from digital tapes of Mariner imagery provided by JPL and compute the mean, 
standard deviation, minimum and maximum of the reflectivities in the array. 

The data required is retrieved from the digital tapes using the Subroutine 
FPIC and subsequently analyzed by the main program whose listing is included 
in Figure A-1. In addition to computing the basic statistics, this main 
program also has the option of computing a histogram of the reflectivity 
values in the specified array. 

Subroutine FPIC retrieves each element in a rectangular array 
located within the digital TV picture data for subsequent processing by the 
main program. The rectangular area of Interest is defined by the following 
calling sequence; 

CALL FPIC (NLS, NSS, NSE, NOL, ARRAY, NLA, NSA, RFAC) where: 

NLS: is the line number within the picture where data retrieval 

begins 

NSS: is the sample number within the line where data retrieval 

begins 
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c 

c 

c 


THIS PROGRAM CALCULATES THE STATISTICS OF A RECTANGULAR IMAGE ARRAY 
SUPPLIED BY SUBROUTINE FPIC. 


ISN 0002 
ISN 0003 
ISN OOOA 
ISN 0005 
ISN 0006 
ISN 0007 


ISN 0008 


ISN 0009 
ISN 0010 
ISN 0011 
ISN 0012 
ISN 0013 


ISN OOIA 


ISN 0015 
ISN 0016 
ISN 0017 
ISN 0018 
ISN 0019 
ISN 0020 
ISN 0021 
ISN 0022 
ISN 0023 
ISN 002A 
ISN 0026 
ISN 0027 
ISN 0029 
ISN 0030 
ISN 0031 
ISN 0032 
ISN 0033 
ISN 0034 


ISN 0036 
ISN 0037 
ISN 0038 
I SN 0039 
ISN 0040 
ISN 0041 
ISN 0042 
ISN 0043 
I SN 0044 
ISN 0046 

ISN 0047 
ISN 0048 
ISN 0050 
ISN 0051 
ISN 0052 
ISN 0053 


DIMENSION R(25,25),IHIST(15) 

REAL*8 TARGIO 
DATA IOPT/0/ 

1 F0RMAT(A8,50X, 414, F4. 0,12) 

9 FORMATllHl) 

10 F0RMAT(6X 'TARGET I DENT I FIC ATI 0N= *,A8,' INITIAL LINE N0. = ',I4, 
*• INITIAL SAMPLE N0.=',I4,' FINAL SAMPLE N0.=',I4,' NO. OF LINE 
♦S=',I4 /// 13X 'PERCENT REFLECTANCE: M INIMUM= ' , F6 .3 , ' MAXIMUM=', 
•F6.3,' AVERAGE = ',F6.3, ' STANDARD DEV I ATION= ' , F8 .5 /// 

♦ 13X 'HISTOGRAM; CLASS I NTERVAL = ' , F6 .3 , // 16X 'CLASS NO. ', 

*15tI3,lX) / 16X 'NO. OF PIXELS ',15(13, IX) / 16X 'TOTAL NO. OF P 

♦IXELS=',I5 /////) 

11 F0RMAT(6X 'TARGET IDENTI FICATION= ',A8,' INITIAL LINE N0.=',I4, 
♦' INITIAL SAMPLE N0.=',I4,' FINAL SAMPLE N0.=',I4,' NO. OF LINE 
♦S=',I4 /// 13X 'PERCENT REFLECTANCE: M INIMUM= ' , F6 .3 , ' MAXIMUM=', 
♦F6.3,' AVERAGE = ',F6.3, • STANDARD DE V I AT ION=' , F8.5 // 13X 'TOTAL 
♦NO. OF PIXELS=',15 /////) 

ICASE=-1 

NLA=25 

NSA=25 

50 READ(5,1,END=100) TARGID ,NL S,N SS ,NSE , NOL , RF AC , lOPT 
ICASE=ICASE^1 
C 

C OBTAIN REFLECTANCE DATA ARRAY 'R' 

C 

CALL FPIC(NLS,NSS,NSE,NOL,R,NLA,NSA,RFAC) 

C 

C COMPUTE THE AVERAGE, STANDARD DEVIATION, MINIMUM AND MAXIMUM OF 

C THE 'R' ARRAY 

C 

SUM=0.0 

SS0=0.0 

RMIN=100. 

RMAX=0.0 
N0S=NSE-NSS>1 
NTOT=NOS*NOL 
DO 60 I>:1,N0L 
DO 60 J=1,N0S 
RT=R( J,I) 

IF(RT.GE.RMIN) GO TO 55 
RMIN=RT 

55 IFIRT.LE.RMAX) GO TO 59 
RMAX=RT 

59 SUM=SUM*RT 

60 SSQ=SSQ*RT**2 
RAVE=SUM/FLOAT(NTOT) 

RSTD=SORT( (SSO-FLOAT(NTOT)*RAVE’»*2)/FLOAT(NTOT-1 ) ) 

IF( IOPT.EO.O)GO TO 80 
C 

C COMPUTE HISTOGRAM 

C 

1 NI=RFAC*(RMAX-RMIN)/15. 

CINT=FL0AT(NI*1 )/RFAC 
DO 65 1*1,15 
65 IHIST(I)=0 
DO 70 1=1, NOL 
DO 70 J=1,N0S 
N=(R( J,I)-RMIN)/CINT +1.0 
70 IHIST(N)=IHIST(N)+1 

IF(M0D(ICASE,4) .EO . 0 ) WR I TE ( 6,9 ) 

WRITE (6, 10) TARGID,NLS,NSS,NSE,NOL,RMIN,RMAX,RAVE,RSTD, CINT, 

*( 1,1=1,15), ( IH I ST( I) ,I=1,15),NT0T 
GO TO 50 

80 IF(M00(ICASE,4) .EQ.O ) WR I TE ( 6, 9 ) 

WRITE (6,11) targid,nls,nss,nse,nol,rmin,rmax,rave,rstd,ntot 

GO TO 50 
100 STOP 
END 

Figure A-1 IMAGE STATISTICS MAIN PROGRAM LISTING 
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SYMBOL 

INTERNAL STATEMENT 

NUMBERS 



I 

0021 

0023 

0038 

0039 

0040 

0042 

0046 0046 0046 0046 0046 0046 

J 

0022 

0023 

0041 

0042 




N 

00A2 

0043 

0043 





R 

0002 

0014 

0023 

0042 




NI 

0036 

0037 






RT 

0023 

0024 

0026 

0027 

0029 

0030 

0031 

MOO 

0044 

0048 






NLA 

0010 

0014 






MLS 

0012 

0014 

0046 

0050 




NOL 

0012 

0014 

0020 

0021 

0040 

0046 

00 50 

NOS 

0019 

0020 

0022 

0041 




NSA 

0011 

0014 






NSE 

0012 

0014 

0019 

0046 

0050 



NSS 

0012 

0014 

0019 

0046 

0050 



SSO 

0016 

0031 

0031 

0033 




SUM 

0015 

0030 

0030 

0032 




CINT 

0037 

0042 

0046 





EPIC 

0014 







lOPT 

0004 

0012 

0034 





NTOT 

0020 

0032 

0033 

0033 

0046 

0050 


RAVE 

0032 

0033 

0046 

0050 




RFAC 

0012 

0014 

0036 

0037 




RMAX 

0018 

0027 

0029 

0036 

0046 

0050 


RMIN 

0017 

0024 

0026 

0036 

0042 

0046 

0050 

RSTD 

0033 

0046 

0050 





SORT 

0033 







FLOAT 

0032 

0033 

0033 

0037 




ICASE 

0009 

0013 

0013 

0044 

0043 



IHIST 

0002 

0039 

0043 

0043 

0046 



TARG 10 

0003 

0012 

0046 

0050 





LABEL 

DEFINED 

REFERENCES 

1 

000 5 

0012 


9 

0006 

0044 

0043 

10 

0007 

0046 


1 1 

000 8 

0050 


50 

0012 

0047 

0051 

55 

0027 

0024 


59 

0030 

0027 


60 

0031 

0021 

0022 

65 

003 9 

0038 


70 

004 3 

0040 

0041 

80 

0048 

0034 


100 

0052 

0012 



LABEL 

AODR 

LABEL 

ADDR 

LABEL 

ADDR 

LABEL 

ADDR 

50 

000DF6 

55 

000ED8 

59 

OOOEE4 

60 

OOOEEE 

65 

OOOFFA 

70 

001054 

80 

001144 

100 

001104 


•STATISTICS* SOURCE STATEMENTS = 52 .PROGRAM SIZE = ■^608 


Figure A-1 (Cont'd) 
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1 • THIS SUBROUTINE WAS CREATED TO PROCESS THE HARS MARINER TV DATA . IT OOCJOLOO 






2 • 

BASICALLY COPIES A TAPE FILE TO A DISK ON THE 

FIRST CALL 

NAD 

00000200 





3 • 

THEN , 

ACCESSES THE DISK IN DIRECT MODE TO FIND 

SQUARES OF 

the 

00000300 





4 « 

PICTURE AS REQUESTED BY THE CALLING PROGRAM 



00000400 





5 

PRINT 

NOGEN 



00000500 

000000 




6 FPIC 

CSECT 




OOOOObOO 




00000 

7 

USING 

«tl5 



00000700 





8 

SAVE 

U4« 12) » ,FPIC_DON_SP ARROW 



00000800 

COOOIA 

18E0 



14 

LR 

14,13 CHAIN THE SAVE AREAS 



00000900 

0000 1C 

41D0 

F030 

00030 

15 

LA 

13, SAVE 



00001000 

000 0 20 

500E 

0006 

00008 

16 

ST 

13, 8( 14) 



00001100 

0000 24 

50E0 

0004 

00004 

17 

ST 

14, 4( 13) 



00001200 

000 0 28 

47F0 

FS8C 

005 8C 

18 

B 

FPICl 



00001300 





20 


DROP 

15 

00001500 



00030 

21 


USING 

SAVE, 13 

OC001600 

0000 30 



22 

SAVE 

DS 

9D 

00001700 




23 

TDCB 

OCB 

ODNAME»TAPE,OSORG>PS,MACRF«<R ),EODADsEOT,SYNAO«TERR, 

*0000 1300 







BUFNQsO 

00001900 




74 

DDCB 

OCB 

DDNAMEsDISK,DSORG>PS,MACRF*UL) ,E0DAD«E00, SYNADsDERR, 

*00002000 







RECFM»F,BUFN0»0,DEVD»0A 

00002100 




125 

ODCBI 

OCB 

ODNAME«OISK,OSORG*DA,MACRF«(RKC),EODAOsEODI,SYNAD30ERRI, 

,*00002200 







RECFM«F,0PTC0»E,LIMCT® 1000,BUFNO»0 

00002300 

0001 BO 



174 

KEY 

DS 

F 

00002400 

000 1 84 



175 

BUFFER 

OS 

CLIOOO 

00002500 



00001 

176 

TRUE 

EQU 

1 

00002600 



00000 

177 

FALSE 

EQU 

0 

00002700 

000 5 6C 

00 


178 

ENDOFHORS DC 

ALMFALSE) 

00002800 

000 5 70 



179 

TEMP 

OS 

0 

00002900 

000 5 78 



180 

NUM8ER0FLINESPERP1CTURE DS F 

00003000 

00057C 



181 

NUHBEROFSAMPLESPERLINE OS F 

00003100 

000580 



182 

nooflabelrecords os f 

00003200 

COO 5 84 



183 

NUMBEROFRECORDSONTAPE DS F 

00003300 

000588 

42000000 


184 

FLOATPAT 

DC 

X *42000000 • 

00003400 


00058C 






186 

FPICl 

OS 

OH 


00003600 

00058C 

1881 





187 


LR 

8,1 

SAVE THE PARM ADDRESS 

00003700 

0005 8E 

9501 

D53C 


0056C 


188 


CLI 

ENDOFHORS, TRUE 

FIRST ENTRY 7 

00003800 

000 5 92 

4780 

D6FE 



0072E 

189 


BE 

FP1C2 

NO > BRANCH 

00003900 







191 


OPEN 

ITOCB, (INPUT)) 

YES - OPEN THE INPUT TAPE 

00004100 






00000 

197 


USING 

IHA0C8,1 


00004200 

0005A2 

4110 

0048 



00078 

198 


LA 

l.TDCS 


00004300 

000 5 A6 

9110 

1030 


00030 


199 


TM 

OCBOFlGS,X* 10* 

WAS OPEN ALRIGHT ? 

000)^4400 

0005 AA 

47E0 

D836 



00866 

200 


BNO 

TOERR 

NO « BRANCH 

00004500 

0005AE 

4140 

OOAO 



00000 

202 


LA 

4 , DOC B 

YES “ GET THE S OF THE DISK OCB 

00004700 







203 


DROP 

1 


00004800 

0005B2 

1B22 





205 


SR 

2,2 

CLEAR THE RECORD COUNTER REGISTER 

COU05000 

0005B4 

47F0 

0596 



005C6 

206 


8 

REAOTAPE 

START THE COPY 

00005100 

0005 B8 






206 

COPYTODISK OS i 

OH 


00005300 







209 


CHECK 

CODECS 

SEE IF DISK WRITE WAS ALRIGHT 

0000 5400 

000 5 C6 






215 

REAOTAPE OS 

OH 


OC005600 

C005C6 

0703 

D5A0 

D5A0 

005D0 

00 5 DO 

216 


XC 

TDECB ,TDECB 

CLEAR THE OECB FOR A READ 

0OUO570J 







217 


READ 

T0ECB,SF,T0CB, 

BUFFER^A.'S* READ A TAPE RECORD 

OuOOSBOO 

0005EE 

4122 

0001 



00001 

230 


LA 

2,112) 

COUNT THE RECORD 

00005900 

000 5 F2 

9501 

D53C 


0056C 


232 


CLI 

ENDOFHORS, TRUE 

ARE WE IN THE DATA RECORDS ? 

00006100 

000 5 F6 

4780 

D624 



00654 

2 33 


BE 

DATARECORO 

YES - BRANCH 

00006200 







235 


CHECK 

TDECB 

NO « SEE IF THE TAPE READ WAS ALRIGHT 

00006400 

000608 

4110 

0158 



00188 

240 


LA 

1,BUFFER*4 

SET BASE FOR USING 

00006500 

00060C 

5920 

09 3 8 



00968 

241 


C 

2 ,=F» 1* 

IS this the first label REOCRD ? 

COU06600 

0006 10 

4720 

0600 



00630 

242 


BH 

NOTISTLABEL 

NO - BRANCH 

00006700 






00000 

244 


USING 

LABEL, 1 

YES - GET NS AND NL FOR THIS TAPE 

OG006900 

000614 

F273 

D540 

1020 

00570 

00020 

245 


PACK 

TEHP.NL 

CONVERT FROM EBCDIC TO BINARY 

00007000 

00061A 

4F00 

05*0 



00570 

246 


cve 

O.TEMP 


00007100 

00061E 

5000 

0548 



00570 

247 


ST 

o.numrsrdflinespfrpicture 

00007200 

000622 

F273 

0540 

1024 

00570 

00024 

249 


PACK 

TEMP , NS 

REPEAT FOR NS 

00007400 

000 6 28 

4F00 

0540 



00570 

250 


CVB 

O.TEMP 


O0UO75O0 

00062C 

5000 

D54C 



0057C 

2 51 


ST 

O.NUMREROF SAMPLE SPER LINE 

00007600 







252 


DROP 

1 


00007700 







254 

* THE 

LAST LOGICAL LABEL RECORD HAS AN *L* IN ThE LAST BYTE - 

00007900 







255 

• 

NO additional label 

RECORDS WILL FOLLOW THIS ONE 

00008000 

000 6 30 






256 

NOTISTLABEL OS 

OH 


OCOOalOO 

0006 30 

4100 

0005 



00005 

257 


LA 

0,5 


00006200 
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0006 34 




2 59 

CHECKF0RLASTLA3EL DS OH 



00008400 




00000 

260 

USING 

LABELRECfl 



0Chj08500 

000634 

9503 

1047 

00047 

261 

CLl 

ENDOFLABELCHECK»C*L 


!S THIS the last LABEL RECORD 

7 00008600 

000638 

4780 

0618 

00648 

2 62 

BE 

SETLASTLABELRECORO 



C0008700 

00063C 

4U1 

004 8 

00046 

264 

LA 

l#72m NO- 

BUNP 

INDEX 

00008900 

000 6 40 

4600 

0604 

00634 

2 65 

8CT 

OtCHECRFORLASTLABEL 

LOOK 

ALLL LABEL RECORDS IN BLOCK 

00009000 

000644 

47F0 

0596 

005C6 

266 

e 

READTAPE GET 

NEXT 

LABEL RECORD 

00U09100 


000648 268 SETLASTLA3ELRECORD DS OH 00C09300 


000648 

9201 

D53C 

0056C 

269 

HVl 

ENDOFHDRS.TRUE yes - SET SWITCH 

00009400 

00064C 

5020 

0550 

00580 

270 

ST 

2»N00FLAeELREC0RDS SAVE THE NUMBER OF LABEL RECORDS 

00009500 

000650 

47F0 

0596 

005C6 

271 

B 

REAOTAPE CONTINUE 

00009600 





272 

DROP 

1 

00009700 


000 6 54 





274 

DATARECORD DS 

OH 


000099U0 






275 



CHECK 

TOcCB 

SEE IF TAPE READ WAS ALRIGHT 

0001*0000 





ooooo 

281 



USING 

IHADC9»4 


00010200 

000 6 62 

9110 

4030 

00030 


282 



TM 

OCBOF LGStXMO* 

IS THE DISK 0C8 OPEN 7 

00010300 

000666 

4710 

0656 


00686 

283 



90 

DDCBO 

YES - BRACNH 

00010400 

00066^ 

5800 

054C 


005 7C 

2 85 



L 

OtNUMBEROFSAMPLESPERLiNE NO > SET THE 8LKSIZE 

00010600 

000 6 6E 

4000 

403E 


0003E 

2R6 



STH 

OtOCBBLKSI 


OOulOTOO 






2sa 



OPEN 

(DOCBtlOUTPUTM OPEN THE OCB 

00010900 

000 67E 

9110 

4030 

00030 


294 



TH 

DCBOFLGS»XMO* 

WAS THE OPEN ALRIGHT 7 

uOOllOOO 

000 682 

47E0 

0846 


00876 

295 



BNO 

OOERR 

NO - BRANCH 

OCOlllOO 

000686 





2 97 

ooceo 

OS 

OH 


00011300 






298 


THE 

KEY OF EACH RECORD IS THE RECORD NUMBER OR LINE NUMBER hHICH 

0001 1400 






299 

• 


RANGES FROM 1 TO NUM8ER0FL INESPERP ICTUR E 

00011500 

000686 

1802 




300 



LR 

0t2 

COMPUTE The KEY 

0001 1600 

000688 

5900 

D550 


00580 

301 



S 

OtNOOFLABELRECOROS 

OC011700 

000 6 SC 

5000 

01 54 


00184 

302 



ST 

OfBUFFER 

SET THE KEY IN FRONT OF THE RECORD 

00011800 

000690 

0703 

D66C 066C 

0069C 

0069C 

304 



XC 

COOECBtCOOECB 

CLEAR THE 0EC6 FOR DISK WRITE 

OOOlWJOtj 






305 



WRITE 

CDOECBtSF.DOCe 

BUFFER WRITE THE RECORD TO DISK 

00012100 

C006BA 

47FF 

068E 


006RE 

319 



8 

*♦41 15) 

CHECK THE RETURN CODE FROM WRITE 

00012300 

0006BE 

47F0 

0586 


00598 

320 



9 

COPYTOOISK 


00012400 

0006C2 

47F0 

0588 


00588 

321 



8 

COPYTOOISK 


00012500 

0006C6 

47F0 

0588 


005B8 

322 



B 

COPYTOOISK 


00012600 






324 



ABEND 

69, DUMP 

WRITE WAS NO GOOD 

00012800 






332 



DROP 

4 


0C012900 

C00 6 0A 





3 34 

EOT 

DS 

OH 


00013100 

000 6 DA 

0620 




335 



BCTR 

2,0 

(2) s NUMBER OF RECORDS UN TAPE 

00013200 

0006 X 

5020 

0554 


00584 

336 



ST 

2,NUHBEROFRECOROSONTAPE 

00U13300 

OO06E0 

SB20 

0550 


00580 

338 



S 

2,NOOFLAR£LRECORDS DOES IT CHECK WITH NL IN HEADER 7 

0CO13500 

0006E4 

5920 

0548 


00578 

3 39 



C 

2 fNJMBEROF L INESPERP ICTURE 

00013600 

0006E8 

4780 

D6CA 


006FA 

340 



BE 

TOOK 

YES • BRANCH 

00013700 






342 



ABEND 

169, DUMP 


00013900 

0006 FA 





351 

TOOK 

DS 

OH 


00014100 






352 



CLOSE 

(TDC9,,0XFI) 

CLOSE THE OCRS 

00014200 

00070A 

4120 

DOE 8 


00118 

36l 



LA 

2,00CB1 

SET DATA TO READ DISK 

00014400 





OOOOO 

362 



USING 

IHA0CB,2 


00014500 

000 70E 

5800 

D54C 


0057C 

364 



L 

0,NUMBER0F SAMPLE SPERLINE 

00014700 

000712 

4000 

2052 


00052 

365 



STH 

0,0CBLRECL 


00014800 

000716 

4000 

203E 


0003E 

366 



STH 

0,0CBBLKS1 


00014900 






368 



OPEN 

(D0CBI,< INPUT)) READY DISK FOR INPUT MODE 

00015100 

000726 

9110 

2030 

00030 


375 



TM 

DC30FLGS,X*10' 

' WAS RE-OPEN ALRIGHT 7 

00015300 

00072A 

47E0 

0856 


00686 

376 



BNO 

DOOERR 

NO - BRANCH 

00015400 






377 



DROP 

2 


00015500 






3 79 

4 

WE 

ENO UP AT 

THIS POINT AFTER THE TAPE IS COPIED TO THE DISK ON THE 

00015700 






380 

• 


FIRST 

CALL, Ok IMMEDIATELY ON THE SECOND AND SUBSEQUENT CALLS 

uOO 15800 

00072E 





381 

FPIC2 

OS 

OH 


00015900 





OOOOO 

382 



USING 

PARMLIST,8 


00C16000 

00072E 

9825 

8000 


OOOOO 

383 



LM 

2»5,NLS 

GET THE ADDRESSES OF THE PARMS 

00016100 

0007 32 

5822 

0000 


OOOOO 

384 



L 

2,012) 


OOOI 02 OO 

0007^6 

5833 

0000 


OOOOO 

385 



L 

3f0(3) 


00016300 

0007 3A 

5844 

0000 


OOOOO 

336 



L 

4,0(4) 


00016400 

000 7 3E 

5855 

0000 


ooooo 

307 



L 

5,0(5) 


00016500 

000 742 

1222 




389 



LTR 

2,2 

IS START LINE NUMBER POSITIVE 7 

00016700 

000 744 

4700 

D386 


008 B6 

390 



BNP 

A269 

NO - BRANCH 

00016800 

C00748 

1812 




392 



LR 

1 ,2 


00017000 

00074A 

1A15 




393 



AR 

1 tS 


00017100 

00074C 

5910 

0548 


00578 

394 



C 

WNUHBEROFLINESPtRPICTURE IS END LINE • IN PICTURE ? 

00017200 

000750 

4720 

0896 


008C6 

395 



BH 

A 369 

NO - BRANCH 

00017300 

000754 

1233 




397 



LTR 

3,3 

IS THE START sample POSITIVE 7 

00017500 

000 7 56 

4 700 

08A6 


00BD6 

398 



BNP 

A 469 

NO - BRANCH 

00017600 

000 7 5A 

1934 




400 



CR 

3,4 

IS START sample less THAN ENO SAMPLE 7 

00017300 

00075C 

4780 

0686 


008 E 6 

401 



9NL 

A 569 

NO - BRANCH 

00017900 

000 7 60 

5940 

054C 


0057C 

403 



C 

4,num3erofsamplesperline is end sample on a line 7 

00018100 

6 

r- 

8 

O 

4720 

09C6 


008F6 

404 



BH 

A 669 

NO - BRANCH 

00018200 


Figure A-2 (Cont'd) 


75 


0007 60 

1244 



406 

LTR 

4,4 

IS END SAMPLE WITHIN A LINE ? 

00016400 

000 7 6A 

4700 

0806 

00906 

407 

8NP 

A 769 

NO - BRANCH 

00018500 

000766 

1256 



409 

LTR 

5,5 

IS the number of lines in square POSITIVE 

00018700 

000770 

4700 

0866 

00916 

4 10 

BNP 

A 869 

NO - BRANCH 

OU018800 

000774 

5860 

8014 

00014 

412 

L 

6, NLA 

GET TH6 DIMENSIONS OF RECEIVING ARRAY 

00019000 

C00778 

5866 

0000 

00000 

413 

L 

6,0(6) 

« OF columns in RECEIVING ARRAY 

OOOlvLOO 

00077C 

5870 

8018 

00016 

414 

L 

7,NSA 


00019200 

000730 

5877 

0000 

00000 

415 

L 

7,0(7) 

« OF ROWS, E.G. ARRAY »A* IS 'A(NSA,NLA)* 

00019300 

000784 

1814 



417 

LR 

lf4 


00019500 

000786 

1813 



418 

SR 

1 ,3 


OOC 19600 

000786 

4111 

0001 

00001 

419 

LA 

1,1(1) 

(1) IS THE NUMBER OF SAMPLES TO PUT IN A 

0001V700 





4 20 

* 


PUT IN A ROW OF >A* ARRAY 

00019600 

00078C 

1917 



421 

CR 

1,7 

IS THE REQUEST FITTaBLE IN ARRAY ROWS 7 

00019900 

000786 

4720 

08P6 

00926 

4 22 

BH 

A 969 

NO - BRANCH 

00020000 

000792 

1956 



4 24 

CR 

5,6 

IS REQUEST FITTABLE IN ARRAY COLUMNS 7 

00020200 

000794 

4720 

0906 

00936 

4 25 

BH 

A 1069 

NO - BRANCH 

00020300 

000 7 98 

5810 

801C 

OOOIC 

427 

L 

1,RFAC 

GET THE SCALE FACTOR 

00020500 

00079C 

7621 

0000 

00000 

428 

LE 

2,0(1) 


00020600 

0007AO 

58C0 

8010 

00010 

430 

L 

12, A 


00020800 

0007A4 

18A7 



431 

LR 

10,7 


C0020900 

0007 A6 

89A0 

0002 

00002 

432 

SLL 

10,2 

NSA*4 - SIZE OF ROW IN ARRAY 

00021000 

000 7 A A 

181A 



4 34 

LR 

1,10 


00021200 

0007AC 

0650 



4 35 

BCTR 

5,0 


00021300 

000 7 A6 

1C05 



436 

HR 

0,5 


00021400 

000 7 80 

4181 

COOO 

00000 

4 37 

LA 

11,0(1,12) 

A ♦ ((NSA-l)*4)*(N0L“l) “ W OF LAST ROW IN 

00021500 





438 

* 


ARRAY TO SET FROM DATA 

00021600 

000784 

0630 



440 

BCTR 

3,0 

READY (3) FOR LOOP AT SAMPLELOOP 

00021800 

000 7 86 

4160 

0001 

00001 

442 

LA 

6,1 

INCREMENT FOR SAMPLE LOOP 

00022000 

000 7 BA 

0640 



443 

BCTR 

4,0 


00022100 

C007BC 

4174 

0154 

00184 

444 

LA 

7,BJFF6R(4) 

a OF LAST SAMPLE IN R60CR0 TO PROCESS 

00022200 

C007CO 




446 

LINELOOP OS 

OH 


00022400 

0007C0 

5020 

0150 

00180 

447 

ST 

2,K=Y 

SET OF RECORD WE WANT TO READ 

00022500 

0007C4 

0703 

07A0 D7A0 

00700 007D0 

449 

XC 

FOECB,FDEC? 

READY FOR FETCH READ 

00022700 





4 50 

READ 

FDEC3,DK,00CBI,3UFP6R,«S', KEY, TEMP ^ETCH DATA FROM DISK 

00022B00 





4 64 

CHECK 

F0EC8 

SEE IF THE READ WAS ALRIGHT 7 

00022900 

000300 

1811 



470 

SR 

1,1 

(1) IS INDEX TO ROW IN A ARRAY 

00023100 

000 8 02 

4143 

0154 

00184 

471 

LA 

4,8UFF6R(3) 

(4) IS PTR TO DAZTA ITEM TO PROCESS 

00023200 

0008 06 




473 

SAMPLFLOOP OS 

OH 


00023400 

000806 

4304 

0000 

00000 

4 74 

IC 

0,0(4) 

GET THE DATA 

00023500 

0008 OA 

4200 

0559 

00589 

4 75 

STC 

0,FLGATPAT^1 

FLOAT IT 

00023600 

000 8 0E 

7800 

0558 

00588 

476 

LE 

o,floatpat 


00023700 

000812 

3002 



477 

DER 

0,2 

SCALE IT 

OC023300 

000814 

7001 

COOO 

00000 

478 

STE 

0,0(1,12) 

STORE IT IN *A* ARRAY 

00023900 

000818 

4111 

0004 

00004 

4 79 

U 

1,4(1) 

BUMP ROW PTR 

00024000 

0008 1C 

8746 

0706 

00806 

480 

BXLE 

4,6,$AMPL6L00P finish THE ROW 

0002**i00 

0008 20 

4122 

0001 

00001 

482 

LA 

2,1(2) 

BUMP DESIRED RECORD KEY 

00024300 

000824 

87CA 

0790 

007CO 

483 

BXLE 

12, 10, LINELOOP process ALL LINES IN SQUARE 

G002440C 


000 8 28 5 6 00 0004 

00004 

485 

L 13,4(13) 

00024600 



486 

RETURN (14,12),RC«0 RETURN 

00024700 


000836 

491 

EOO 

DS 

OH 

00024900 


492 


ABEND 

1169, DUMP 

00025000 

C00 8 46 

501 

DERR 

OS 

OH 

0002520U 


502 


ABEND 

1269, DUMP 

00025300 

000856 

511 

TERR 

DS 

OH 

00025500 


512 


ABEND 

1369, DUMP 

00025600 

000 8 66 

521 

TOERR 

OS 

OH 

00025800 


522 


ABEND 

1469, DUMP 

000^5900 

LOC OBJECT CODE 

AOORl A00R2 stmt 

SOURCE 

STATEMENT 

ASM H V 05 11.17 07/14/77 

000 8 76 

531 

OOERR 

DS 

OH 

00026100 


532 


ABEND 

1569, DUMP 

00026200 

000 8 86 

541 

OOOERR 

DS 

OH 

00026400 


542 


ABEND 

1669, DUMP 

00026500 

0008 96 

551 

DERRI 1 

OS 

OH 

00026700 


5 52 


ABEND 

1769, DUMP 

00026300 

C003 A6 

561 

0ERRI2 

DS 

OH 

00027000 


562 


ABEND 

1869, DUMP 

00027100 
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nooaP6 

571 A269 

OS 

OH 

00027300 


5 72 

ABEND 

269, DUMP 

0OC274O0 
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581 A369 

OS 

Oh 
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000 8 06 
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0O0279OU 


592 

ABEND 

469, DUMP 

00028000 

0008E6 

601 A569 

OS 

OM 

0C02b20O 


602 

ABEND 

569, DUMP 

00023300 

0008F6 

611 A669 

OS 

OH 

0002B500 


612 

ABEND 

669, DUMP 
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0269 
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NSE; is the last sample within the line to be included in the 


retrieved data 

NOL: is the number of consecutive lines for retrieval of data 

(Including NLS) 

These parameters define the rectangular area within the image. 
ARRAY: is the array into which the retrieval data is put 

NLA: is the second dimension of the array (the number of 

rows (lines) the array contains) 

NSA: is the first dimension of the array (the number of columns 

(samples) the array contains) 

These parameters define the array to receive the data. 

RFAC; is a scale factor by which each san^le is divided 
Note; NLS, NSS, NSE, NOL, NLA, NSA are assumed to be INTEGER *4 
arguments while ARRAY and RFAC are assumed to REAL *4 
arguments and further, a previous statement: 

DIMENSION ARRAY (NSA, NLA) is assumed to have been included in 
the program that uses FPIC, which is FORTRAN compatible. 

The data to be manipulated by FPIC is assumed to reside on a tape 
defined by a DD card as follows; 

//GO. TAPE DD DSN=name,DISP=OLD,UNIT=2400, 

// V0L=SER=vol ser ,DCB= (RECFM=U ,BLKSIZE=1000) 

// LABEL=(flle #,SL,,IN) 

As part of FPIC's operation, another DD card is required as follows 
//GO .DISK DD DSN=&&TEMP,DISP=NEW,UNIT=SYSDA, 

// DCB =(DSgTRG=DA,KEYLEN=4, SPACE =(TRK, (300,10) ) 
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The operation of FPIC is as follows: 

1) The first call to the subroutine causes the data (not 
the label) records to be copied from the tape to the disk 
data set as a direct access file with the line number as 
the record key. 

2) The required lines (records) are then read from the disk, 
the sample bytes are extracted, floated, scaled and stored 
in ARRAY. 

Certain operational errors will cause user abends. These are listed below. 


Abend 

Reason 

69 

disk write rejected 

169 

no. of data records NL in header 

269 

start line number negative 

369 

requested number of lines >NL in header 

469 

start sample # negative 

569 

start sample # > end sample # 

669 

end sample # ^ NS in header 

769 

end sample # negative 

869 

number of lines in square (NOL) negative 

969 

(end sample # - start sample #)>NS in A array 

1069 

number of lines (N0L)> NL in A array 

1169 

EODAD on disk (write) 

1269 

disk I/O error (write) 

1369 

Tape I/O error 

1469 

unable to open tape data set deb 
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1569 


unable to open disk data set (write) deb (write) 
1669 unable to open disk data set deb (read) 

1969 end of dataon disk data set (read) 

2069 I/O error on disk data set (read) 

A program was written to eompute the average speetral refleetanee 
within a partleular speetral band given the speetral refleetlvlty of the 
sample and the speetral transmission of the filter under e onsi derat ion . 

Figure A-3 is a listing of this program whleh was written in Fortran language. 
Eaeh of the steps in the program is elearly Identified by eominents and no 
further explanation is required. As utilized in the present study, FIL (1,J) 
through FIL (6,J) represent the six speetral transmissions for the blue, green, 
red filters of Mariner VI and Marlnter VII r espeetively . REAR represent the 
average speetral refleetlvlty of the sample in eaeh of these speetral bands 
and the value of RATIO (1) and RATIO (2) are the red-to-blue and red-to-green 
refleetlvlty ratios for the Mariner VI speetral bands while RATIO (3) and 
RATIO (4) are the eorrespondlng values for the Mariner VII speetral bands. 


81 


I SN 0002 
ISN 0003 
ISN 000^ 
ISN 0005 


ISN 0006 


ISN 0007 


ISN 0008 
ISN 0009 
ISN 0010 
ISN 0011 
ISN 0012 
ISN 0013 


ISN 0014 
ISN 0015 
ISN 0016 
ISN 0017 
ISN 0018 
ISN 0019 
ISN 0020 


ISN 0021 
ISN 0022 
ISN 0023 


C 

c 

c 

c 

c 

c 

c 


:^t************m***m*************** ********************* *********************** 
** ** 

♦♦ THIS PROGRAM COMPUTES THE AVERAGE SPECTRAL REFLECTANCE OF A SAMPLE 
** FOR SPECTRAL BANOS GIVEN THE SPECTRAL TRANSMISSION AND THE *♦ 

** SAMPLE REFLECTANCE. 

** ** 
***********m******* ************************************ *********************** 

DIMENSION FIL(6,151 ) ,R( 1 51 ) ,S I D( 20 ) ,RB AR ( 6) .RAT 10 <4 ) 

1 F0RMAT(4( lOX.FlO.O) ) 

2 FORMAT(20A4) 

3 FORMAT(1H1,9X,20A4//15X 'AVERAGE REFLECTANCE*/ 20X • B AND' , 2X , • 1 • . 

■*8X, '2',8X,'3',8X,'4',8X, '5',8X,'6'// 20X, 6F9.3 //15X 'REFLECTANCE 

* RATIOS', lOX, 4F8.3, 4(/)) 

4 FORMAT! lOX , 20A4//1 5X 'AVERAGE REFLECTANCE*/ 20X ' B AND ' , 2X, ' 1 • , 

* 8X,'2',8X,'3*,8X,'4',8X, '5',8X,'6'// 20X, 6F9.3 //15X 'REFLECTANCE 

* RATIOS', lOX, 4F8.3, 4(/l) 

CALL CLEAR(FIL( 1,1),FIL(6,151) ) 

****m************* ******************************************** *************** 
* * 
* READ SPECTRAL TRANSMISSION DATA * 

« * 
****************** *********************************************************** 


C 

c 

c 

c 

c 

c 

c 


READ(5,1) 
READ! 5,1) 
READ(5,1 ) 
REA0(5,1) 
READ!5,1) 
READI 5,1) 


!FILI1,J) ,J=2,151) 
(FIL(2,J) ,J=2,151) 
(FIL(3,J) ,J=2,151) 
(FIL(4,J) ,J=1,151) 
(FIL(5,J) ,J=1 ,151) 
(FIL(6,J),J=1,151) 


C 

c 

c 

c 

c 

c 

c 


***************************************************************************** 
* * 

* NORMALIZE SPECTRAL TRANSMISSION DATA ♦ 

* * 
***************************************************************************** 

DO 30 1=1,6 
S=0.0 

DO 22 J=l,151 
22 S = S + FIUI,J) 

DO 24 J=l,151 
24 FIL(I,J)=FIL(I,J)/S 
30 CONTINUE 

***************************************************************************** 
* * 

* READ SAMPLE SPECTRAL REFLECTANCE DATA ♦ 

* * 
***************************************************************************** 

ICASE=0 

40 READ(5,2,END=70) SID 
iICASE=ICASE+l 


Figure A-3 AVERAGE REFLECTIVITY MAIN PROGRAM LISTING 
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ISN 0024 


REAO(5,l) R 


ISN 0025 
ISN 0026 
ISN 0027 
ISN 0028 
ISN 0029 


I SN 0030 
ISN 0031 
ISN 0032 
ISN 0034 
ISN 0036 
ISN 0038 


ISN 0040 
ISN 0042 
ISN 0043 
I SN 0044 
ISN 0045 


ISN 0046 
ISN 0047 


C 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


* * 

* DETERMINE AVERAGE SAMPLE REFLECTANCE FOR EACH SPECTRAL BAND ♦ 

* * 


DO 50 1=1,6 
RBAR( n=0.0 
DO 45 J=l,151 

45 RBAR( I)=RBAR( I)+R(J)*FIL(I,J» 
50 CONTINUE 


4 m 

* COMPUTE RATIOS OF AVERAGE REFLECTANCES ♦ 

* * 
m**************************************************************************** 


DO 55 N=l,4 
55 RATI0(N)=0.0 

IF<RBAR( 1 ) .GT.O.OJ R ATIO ( 1) =RB AR ( 3 )/RB AR( 1) 

IF(RBAR(2).GT.0.0) R ATI 0 ( 2 > =RB AR ( 3 )/RB AR( 2) 

IF»RBAR(4).GT.0.0) R ATIO < 3 ) =RBAR ( 6 ) /RB AR ( 4> 

IF(R3AR(5 ).GT.O.O) R AT 10 (4 ) =RBAR ( 6 )/RB AR ( 5) 

m^f^m********m**************m****************** ******************************* 
* * 

* PRINT OUTPUT DATA AND RECYCLE TO PROCESS NEXT CASE, IF ANY. * 

* * 
***************************************************************************** 


IF(M0D(ICASE,5) .NE.l ) GO TO 60 

WRITE ( 6 , 3) (SID(K) ,K = 1,20) , (RBAR( I ) ,I = 1,6 ) , ( RAT 10 IN ) , N= 1,4 ) 

GO TO 40 

60 WRITE(6,4) I SID (K ) ,K = 1 , 20) , I RB AR 1 1 ) , I = 1 ,6 ) , ( RAT 10 (N ) , N= 1 , 4 ) 

GO TO 40 

***************************************************************************** 
* * 

♦ TERMINATE EXECUTION ♦ 

• * 
*************************************************************************: •f*tr 

70 STOP 
END 
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Figure A-3 (Cont'd) 
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